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Abstract 

This  paper  introduces  and  analyzes  an  efficient  algorithm  for  the  random 
accessing  of  a  broadcast  channel  by  a  large  number  of  packet-transmitting, 
bursty  users.  The  algorithm  uses  a  mini-slot  overhead,  per  packet,  to  extract 
detailed  information  regarding  possible  packet  collisions.  In  the  event  of  a 
collision,  this  information  is  used  by  the  algorithm  for  accelerating  the 
collision  resolution  process.  The  maximum  stable  throughput  and  the  mean 
packet  delay  induced  by  the  algorithm  are  evaluated  via  a  systematic  analysis 
method.  The  packet  delay  characteristics  indicate  that  the  proposed  algorithm 


is  well  suited  for .operation  on  satellite  channels,  over  a  wide  range  of  input 
traffic  rates.  , 
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1.  INTRODUCTION 


We  consider  the  random-accessing  of  a  broadcast  communication  channel, 

(such  as  a  satellite  link,  or  a  ground  radio  channel) ,  by  a  large  number  of 
independent,  packet-transmitting,  bursty  users,  whose  cumulative  'jacket 
generating  process  is  Poisson.  Overlapping  transmissions  result  in  a  collision, 
and  all  the  packets  involved  must  be  retransmitted.  Users  monitor  the  channel 
activity  and  acquire  some  form  of  feedback  information  associated  with  packet 
transmissions.  The  feedback  information,  whose  specific  form  depends  on  the 
precise  model  at  hand,  is  then  used  by  the  users  for  scheduling  packet  transmissions 
and  retransmissions  in  accordance  with  a  distributed  control  scheme,  referred 
to  as  random-access  algorithm  (RAA) . 

Given  a  specific  model  for  the  available  feedback  information.,  the  performance 
of  a  random-access  algorithm  is  commonly  measured  in  terms  of  the  packet  delay 
characteristics  that  it  induces,  and  the  maximum  throughput  that  it  attains. 

For  the  Poisson  infinite-population  user  model,  the  maximum  throughput  is  defined 
as  the  maximum  input  rate  that  the  algorithm  maintains  with  finite  mean  packet 
delay.  A  random-access  algorithm  is  called  stable  if  its  maximum  throughput  is 
greater  than  zero. 

The  design  and  performance  of  a  random-access  algorithm  depend,  critically, 
on  the  availability,  quality,  and  timely  acquisition  of  the  feedback  information, 
which  in  turn  depend  on  the  channel  and  user  characteristics  of  the  specific 
application  at  hand.  Among  the  several  aspects  that  characterize  the  feedback 
information  structure  we  distinguish  the  following:  1)  feedback  level,  2)  feedback 
delay.  The  feedback  level  specifies  what  kind  of  information,  regarding  the 
transmission  activity  on  the  channel,  can  be  extracted  by  the  users  who  monitor 
the  channel.  The  feedback  delay  specifies  when  this  information  becomes  available 
to  the  users,  and  is  defined  as  the  time  it  takes  for  a  user  to  determine  the 
current  transmission  activity  on  the  channel.  To  be  specific,  let  us  assume  that 


packets  are  of  equal  length,  the  channel  is  slotted  with  slot  size  equal  to  the 
packet  transmission  time,  and  that  the  users  initiate  transmission  only  at  slot 
boundaries.  In  this  case,  the  feedback  information  is  commonly  modeled  as  providing 
the  "outcome"  of  each  slot,  at  the  end  of  each  slot,  (i.e.,  the  feedback  delay  is 
assumed  equal  to  one  slot),  where  various  such  outcomes  may  be  distinguished, 
depending  on  the  feedback  level.  The  most  thoroughly  studied  feedback  level  is 
the  ternary,  which  distinguishes  among  empty,  busy  with  one  packet  (or  successful), 
and  busy  with  at  least  two  packets  (collision)  slots.  It  has  been  shown  that,  under 
ternary  feedback,  the  achievable  maximum  throughput  cannot  exceed  0.587  packets  per 
slot  [1];  the  most  efficient  algorithm  known  to  date  attains  a  maximum  throughput 
of  0.487  packets  per  slot  [2].  In  some  applications,  due  to  special  modulation, 
coding  and  encryption  techniques,  or  to  channel  noise  conditions,  ternary  feedback 
might  not  be  easily  available.  In  such  applications  a  binary  feedback  model  is 
more  realistic,  and  includes  the  following  cases  [3]:  Distinguishing  between  either 
collision  versus  noncollision  slots,  (C-NC  binary  feedback)  or  empty  versus  nonempty 
slots,  (E-NE  binary  feedback)  or  successful  versus  nonsuccessful  slots,  (S-NS 
binary  feedback).  Multiary  feedback  models  have  also  been  considered  [4],  where 
it  is  assumed  that  the  existence  of  energy  detectors  enables  users  to  determine 
the  number  of  collided  packets  within  each  collision  slot,  whenever  this  number 
is  below  a  certain  limit;  however,  even  with  this  more  informative  level  of  feedback, 
the  maximum  throughput  of  the  existing  algorithms  does  not  exceed  0.532  packets  per  slot. 

Given  a  feedback  level,  the  performance  of  RAAs  can  be  substantially 
improved  if  the  feedback  delay  is  small  compared  to  the  packet  transmission  time. 

This  is,  for  example,  the  case  in  the  multiple  access  channels  used  in  some  local 
area  networks,  where  simple  RAA's  based  on  carrier  sensing  [5,6],  attain  maximum 
throughputs  close  to  unity  with  uniformly  low  packet  delays. 

For  channels  with  feedback  delay  equal  to  or  larger  than  the  packet  transmiss¬ 
ion  time,  improvements  in  the  maximum  throughput  are  hard  to  obtain,  as  long  as 
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the  algorithms  rely  on  just  the  passive  observation  of  channel  history.  For  such 
channels,  the  maximum  throughput  can  be  improved  by  the  use  of  reservation-based 
algorithms  [7,8].  Reservation  algorithms  restrict  the  multiple  access  contention 
problem  to  a  reservation  subchannel.  The  reservation  subchannel  uses  a  fraction 
of  the  total  channel  capacity  for  the  transmission  of  explicit  reservation  requests 
(mini-packets) ,  and  packets  are  transmitted  collision-free  following  successful 
reservation  requests.  If  the  reservation  overhead  per  packet  is  sufficiently 
small,  reservation  algorithms  can  attain  maximum  throughput  close  to  unity, 
independently  of  how  large  the  feedback  delay  is.  In  channels  where  the  ratio 
between  feedback  delay  and  packet  transmission  time  is  large,  however,  this 
throughput  improvement  is  obtained  at  the  expense  of  higher  mean  packet  delay  at 
low  input  rates,  as  compared  to  the  mean  packet  delay  induced  by  stable  RAA's. 

The  reason  for  this  delay  increase  is  that  packet  transmissions  are  Invariably 
associated  with  a  minimum  delay  equal  to  the  feedback  delay  of  a  reservation  request. 
Considering,  for  example,  a  satellite  channel  where  the  round-trip  propagation 
delay  is  equal  to  270  msec,  the  use  of  reservations  implies  that  packets  suffer 
a  delay  of  at  least  540  msec  before  they  reach  their  destination,  even  for  near 
zero  input  rates.  Note  that  the  additional  delay  Introduced  by  the  required 
transit  time  of  a  reservation  request  may  prove  Intolerable  in  application  involving 
time  critical  interactive  traffic,  such  as  a  remote  computing,  data  base  query,  etc. 

Hybrid  schemes  that  combine  the  better  features  of  random  access  and  reservation 
access  have  also  been  proposed  [9,10].  These  schemes  use  TDMA.  or  some  other  fixed 
assignment  technique  for  the  transmission  of  reservation  packets,  and  for  small 
user  populations,  they  have  excellent  throughput-delay  characteristics.  However, 
their  performance  deteriorates  rapidly  as  the  number  of  users  increases,  since  the 
reservation  overhead  is  proportional  to  the  user  population  size. 

In  this  paper,  we  present  and  analyze  a  RAA  that  utilizes  more  detailed 
feedback  information  than  the  commonly  used  ternary  feedback  to  accelerate  the 
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collision  resolution  process,  and,  therefore,  to  attain  substantially  higher  maximum 
throughput.  As  explained  below,  the  more  detailed  feedback  information  is  created 
usii  •  a  aartial  reservation  strategy,  which  is  based  on  the  exchange  of  imprecise 
resehai  on  requests  conveying  only  a  few  "bits"  of  information;  these  partial 
reservation  requests  require  little  bandwidth  on  the  channel,  and  their  size  is 
independent  of  the  number  of  the  channel  users.  In  contrast  to  reservation 
algorithms,  for  channels  with  long  propagation  delay,  the  proposed  algorithm 
exhibits  good  delay  performance  in  both  the  low-throughput  and  the  mid-throughput 
range. 

To  create  the  more  detailed  feedback  the  algorithm  divides  time  into  contiguous 
slots  of  equal  length,  taken  as  the  unit  of  time.  Each  slot  consists  of  a  control 
slot  (CS)  and  a  data  slot  (DS) .  The  CS  is  divided  into  m  control  mini-slots  (CMS) , 
each  of  which  is  of  one  bit  (or  a  few  bits)  duration;  the  DS  is  long  enough  to 
accomodate  a  standard  data  packet,  (see  figure  1).  Each  time  a  user  decides  to 
transmit  a  packet,  it  transmits  it  in  the  DS  and  simultaneously  it  Injects  a  pulse 
in  one  of  the  m  CMS's  chosen  with  equal  probability.  It  is  assumed  that  at  the  end 
of  slot  i+P  all  users  can  determine  without  error  whether  the  transmissions  in  the 
i-th  slot  resulted  in  a  packet  collision  in  the  DS  (i.e.,  C-NC  binary  feedback) , 
and  whether  a  CMS  of  the  i-th  slot  is  empty  (i.e.,  E-NE  binary  feedback) ;  P  denotes 
the  channel's  propagation  delay  in  slots.  The  additional  feedback  information  provided 
by  the  CS  is  used  by  the  users  to  prevent  idle  slots  and  further  collisions  during 
the  collision  resolution  phase,  and  is  the  principal  mechanism  for  improving  perform¬ 
ance  .  Since  the  feedback  from  each  CMS  provides  information  only  on  the  presence 
of  absence  of  a  pulse,  the  overhead  needed  for  the  GS  is  small.  We  shall  see  that 
this  small  investment  of  channel  capacity  results  in  substantial  Improvement  in  the 
overall  system  performance. 

The  control  mini-slot  concept,  in  conjuction  with  collision  resolution 
algorithms,  has  been  used  by  Huang  and  Xu  [11] ,  where  ternary  feedback  is  assumed 
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available  from  each  mini-slot,  and  Merakos  and  Kazakos  [12],  where  the  mini-slot 
feedback  is  of  the  binary  type  considered  here.  The  algorithm  considered  in  this 
paper  is  a  generalized  version  of  the  second  algorithm  in  [12], and  is  very  similar 
to  the  one  proposed  and  analyzed  recently  by  Huang  and  Berger  [13],  [14],  after  an  in 
dependent  and  concurrent  investigation.  An  important  contribution  of  the  present 


paper  is  a  novel  delay  analysis  method,  which  is  based  on  the  asymptotic  properties 
of  regenerative  processes.  The  proposed  analysis  method  is  of  independent  interest, 
and  can  be  used  for  evaluating  the  performance  of  a  large  class  of  RAA's. 
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The  ZPDA  is  defined  by  the  following  two  rules: 

Rule  1  (First  Tim?  Transmission) 

All  users  observe  the  feedback  at  the  end  of  each  slot  and  utilize  a 
counter,  called  global  counter,  to  determine  when  to  transmit  a  packet  for  the 
first  time.  A  user  initializes  and  updates  his  global  counter  according  to  the 
following  rule. 

The  counter  is  set  to  1  at  the  beginning  of  the  first  slot,  then  incremented 
by  (j-1)  at  the  end  of  each  collision  slot  with  j  nonempty  mini-slots, 
and  decremented  by  one  at  the  beginning  of  each  collision-free  slot.  When 
the  counter  reaches  0  it  is  immediately  reset  to  1  and  the  above  described 
operation  is  repeated  indefinitely. 

Let  denote  global  counter  indication  at  the  beginning  of  slot  i,  i>0, 
and  define  the  sequence  {l^}n>0  on  the  channel  axis,  and  the  sequence  (tn)n>0 
on  the  arrival  axis  as  follows: 

In~l,  I  .  =min{i  :  G.=0,  i>I  }  ,  n=0,l,2,... 
u  n+i  l  n 

V0’  tn+l"tn  +  mln  ‘W4’ 

where  A  is  a  positive  algorithmic  parameter. 

Using  the  above  sequences,  the  first  time  transmission  rule  can  be  expressed 


as  follows : 
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A  packet  arrived  at  time  instant  t  is  transmitted  for  the  first  time  in 

slot  I  if  te  [ t  , t  , ) . 
n  n  n+1 

Rule  2  (Retransmission) 

Each  user  with  a  packet  involved  in  a  collision,  determines  when  to  re¬ 
transmit  his  packet  using  a  second  counter,  called  local  counter,  in  accordance 
with  the  following  rule. 

Following  a  collision  in  which  he  is  involved  the  user  sets  his  local 
counter  to  (i-1) ,  if  he  injected  his  pulse  in  the  i-th  nonempty  mini¬ 
slot.  Then,  he  increments  it  by  (j-1)  for  each  subsequent  collision  slot 
with  j  nonempty  mini-slots,  and  decrements  it  by  one  for  each  subsequent 
collision-free  slot.  When  the  local  counter  reaches  C.  the  user  retransmits 
in  the  next  slot. 

Remarks 

1)  As  it  can  be  seen  from  Rule  1,  the  ZPDA  divides  the  arrival  axis  into 

contiguous  time  intervals  [t  ,t  , ),  n>0.  The  packets  (if  any)  arrived  in  the 

n1  n+i  — 

interval  [t  ,t  are  first  transmitted  in  slot  I  .  If  there  is  a  collision  in 
n  n+1  n 

I  ,  the  collision  is  resolved,  (i.e.,  the  collided  packets  are  eventually 
n 

successfully  transmitted),  during  slots  1^+1,  In+2 , . . . ,In+^-l,  in  accordance  to 
Rule  2.  This  is  illustrated  in  figure  2,  where  a  collision  of  three  packets 
is  resolved. 

The  interval  [tn,tn+^)  will  be  referred  to  as  the  n-th  enabled  arrival  interval 

(EAI);  the  interval  [I  ,1  ..)  will  be  referred  to  as  the  n-th  collision  resolution 

n  n+1  - 

interval  (CRI) ,  and  the  instants  I  ,n>0,  will  be  referred  to  as  collision  resolution 
-  n  —  - 

instants. 

2)  Rule  2  defines  a  collision  resolution  algorithm  (CRA) ,  which  is  a  simple 
modification  of  the  well-known  Tree  CRA  [  18  ] .  The  modification  consists  in 
utilizing  the  feedback  information  provided  by  the  m  CMS's  to  skip  the  empty 
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branches  in  an  m-ary  tree  search,  and,  therefore,  accelerate  the  collision 
resolution  process. 

3)  The  algorithm  utilizes  an  elementary  window  flow  control  mechanism  on  the 
arrival  axis  by  not  allowing  the  length  of  an  EAI  to  exceed  A.  The  window  size 
A  will  be  optimized  for  throughput  maximization. 

2.2  Delay  Analysis  of  the  ZPDA 

Consider  the  random  access  system  operating  with  the  ZPDA  over  the  time 
interval  [0,  -H»).  Let  the  arriving  packets  be  labelled  n=l,2,...,  according 
to  the  order  of  arrival.  We  define  the  delay,  P^,  experienced  by  the  nth  packet 
as  the  time  from  its  arrival  at  the  transmitter  until  the  completion  of  its 
successful  reception  by  the  receiver,  (so  that  P =P+1,  when  the  nth  packet  is 
successfully  transmitted  beginning  at  the  same  moment  it  arrives  at  the  trans¬ 
mitter  -  for  the  ZPDA,  P=0) .  We  are  interested  in  evaluating  the  mean  packet 
delay  induced  by  the  ZPDA. 

From  the  rules  of  the  algorithm  it  can  be  readily  seen  that,  at  each  collision 

resolution  instant,  all  packets  arrived  before  time  t^  have  been  successfully 

transmitted,  and  the.  packets  arrived  in  the  interval  [t  ,1  )have  not  accessed 

n  n 

the  channel  yet.  Let  d  ^  I  -t  ,n>0;  d  will  be  referred  to  as  the  lag  at  I  . 

n  n  n  —  n  — “•  n 

Since  t^O,  Iq*1»  we  have  dg=l.  Let  tq^Iq  and  define  as  the  first  collision 

resolution  instant  after  T . ,  i-0,1,2,...,  at  which  d_  -1. 

1  Ti+1 

From  the  operation  of  the  algorithm  it  can  be  seen,  after  a  little  thought, 
that  the  delay  process  "probabilistically  restarts  itself"  at  each  instant 

T^,  i>0.  The  interval  l^,T^+^),  i>0,  will  be  referred  to  as  the  ith  session; 
since  the  packet  arrival  process  (Poisson)  is  memoryless  and  has  independent 
increments,  the  session  lengths,  T^+j-T^,i>0,  are  independent,  identically 
distributed  (i.i.d.)  random  variables. 


Let  R^,  i>0,  denote  the  number  of  packets  successfully  transmitted  in  the 


interval  (T  ,T  ] .  (Note  that  R^  denotes  also  the  number  of  packets  arrived 


during  (0,1\-1)).  Accordingly,  C^=  R^+^-R^,  t^ie  num^er  °f  packets 


successfully  transmitted  during  the  ith  session.  The  sequence  {R^K>q  is  a 
renewal  process,  since  {CP^q  is  a  sequence  of  nonnegative  i.i.d.  random 


variables.  Furthermore,  the  delay  process  {0  }rt>^  regenerative  with  respect 


to  the  renewal  process  w^th  regeneration  cycle  C^;  that  is,  the  process 


{Or  +n)n>i>  f°r  every  i>0,  is  a  probabilistic  replica  of  the  process  {0n)n>^. 


The  regenerative  probabilistic  structure  of  the  delay  process  makes  it 
possible  to  express  its  asymptotic  behavior  in  terms  of  quantities  that  refer 
only  to  one  regeneration  cycle  of  the  process.  This  is  made  precise  by  the  follow¬ 
ing  elegant  and  powerful  result  from  the  theory  of  regenerative  processes  t  15  1  - 
Theorem  1 


C- 


If  C=E(Cq)<°°  and  S=  E(  then  there  exists  a  real  number  D  such  that 

n 


D=lim  —  .1  V.  =lim  —  E{  E.O. }  =  w.p.l 

n4oo  n  i=l  i  n4*0  n  i=l  i  C 


Furthermore,  if,  in  addition  to  the  finiteness  of  C  and  S,  the  distribution 


of  Cq  is  not  periodic,  then  converges  in  distribution  to  a  random 


variable  V  ,  and 


E  (P  )  =  f 

CO  c 


In  our  case,  the  distribution  of  Cq  is  clearly  not  periodic;  thus,  provided 
that  both  C  and  S  are  finite,  the  limiting  average,  the  limiting  expected 
average,  and  the  mean  of  the  limiting  distribution  of  exist, 

coincide,  and  are  finite;  their  common  value  D  will  be  referred  to  as  the  mean 
packet  delay.  Note  that  D  can  be  computed  using  only  the  per  cycle  quantities 
C  and  S. 

Next  we  develop  two  systems  of  equations,  whose  solution  is  then  used  to 
compute  the  mean  cycle  length  C,  and  the  mean  cumulative  delay  over  a  cycle  S  . 


Consider  first  the  mean  cycle  length  C.  If  the  mean  session  length 


H  =  E{T.+^  -  T.},  i>0,  is  finite,  then  by  Wald's  theorem  we  have  that 

C  =  X  H  (1) 

To  determine  H  we  proceed  as  follows. 

Let  1(d)  be  a  collision  resolution  instant  at  which  the  lag  is  equal  to 
d;  let  T(l)  be  the  first  collision  resolution  instant  after  1(d)  at  which  the 
lag  is  equal  to  one.  The  interval  f  I C<i)  >  T(l))  is  called  a  d-session.  Next 
we  define  the  following  random  variables: 

h,  :  the  length  of  a  d-session 
d 

t,  :  the  length  of  the  first  CRI  of  a  d-session 
d 

Note  that,  by  definition,  H  =  E(h^) , 

From  the  description  of  the  algorithm  we  have  that  the  lags  at  two  successive 
collision  resolution  instants  satisfy  the  following  relation. 


d.  if  d  <  A 

l  i  — 

d  -  A  +1,  if  d  >  A 

i  d .  1 

l 


From  (2)  and  since  d^  =1,  we  deduce  that  d^  e  A,  for  every  i>0t  where  A  is  a 
denumerable  subset  of  the  set  of  positive  real  numbers;  specifically 


A  ■  {d  :  d“k-iA>l;  k,i  integers) 


From  the  operation  of  the  algorithm  and  in  view -of  (2),  it  can  be  readily 
deduced  that  the  h^,  d  e  A,  satisfies  the  following  relations: 


lA  if  l .  =  1 

d  d 

li  *  \  1£  V1 

d 


,l<d<A 


(3.  a) 
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»d  -*d  +  hd-A+£,  ’  d>4  <3'b) 

d 

If  we  let  H,  ■  E  (h .) ,  then  taking  expectations  In  (3)  yields 
d  d 

Hd  -  E  (£d)  +  £?2  pd(l)  Hd  ,  l<d  <  A  (4. a) 

Hd  ‘  E  «4>  +  i?l  Pd(1)  “d-A+l  >.  d  >  4  <4-b> 


where  pd(i)  is  the  probability  distribution  of  £d»deA.  Note  that  for  d  >  A  , 

E  (£d)  «  E  (Ip) ,  and  pd(i)  ■  p^(i) ,  since  when  the  lag  is  greater  than  A  the  enabled 


interval  is  of  length  A.  Thus,  the  sequence  satisfies  the  following  infinite 


dimensional  system  of  linear  equations 


bd  + 


Cdt  *t 


d  e  A 


where  bd  *  E  (t^) ,  l<d<  A,  bd  -  E  (t^ ,  d  >  A  .  and  where  cdfc  ,  d,t  e  A  ,  are 
nonnegative  coefficients  that  can  be  appropriately  identified  from  (4). 

Since  H  ■  H^,- we  have  from  (1)  that  C  ■  X  H^.  Thus,  to  compute  C  we  need 
to  compute  H^.  This  will  be  done  by  solving  system  (5).  Before  examining  the  solu¬ 
tion  of  system  (5),  however,  we  first  develop  a  similar  system  of  equations  that 
can  be  used  to  compute  the  mean  cumulative  delay  S.  For  the  development  of  such 
a  system  we  need  the  following  definitions: 

wd  :  the  sum  of  the  delays  experienced  by  all  the  packets  that  were  successfully 
transmitted  during  a  d-session. 

(t) .  :  the  sum  of  the  delays  experienced  by  all  the  packets  that  were  successfully 

O 

transmitted  during  the  first  CRI  of  a  d-session. 

The  operation  of  the  algorithm  yields  the  following  relations  for  the 
wd>  d  e  A. 
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if  l,  =  1 

a 


w ,  +  u,  if  t,  >1 

a  JL  a 

a 


,  1  <d<  A 


Wd  =  “d  +  Wd-A+£ 


,  d  >  A 


(6.a) 


(6.b) 


If  we  let  =  E  (w^) ,  then  taking  expectations  in  (6)  yields 


Wd  =  E  (o)d)  +  X2  Pd(i)  V±  ,  1  <  d  <  A 


(7. a) 


Wd  E  (“d*  +  i=l  Pd(i)  Wd-A+i  ’ 


d  >  A 


(7.b) 


From  (7)  we  have  that  the  sequence  (W.},d  e  A,  satisfies  the  following  infinite 

d 

dimensional  system  of  linear  equations 

*d  '  bd  +  tiA  cdt  \  •  d  s  A  <8) 

* 

where  b,  =  E  (w,)  ,  d  e  A  ,  and  the  nonnegative  coefficients  c.. ,  d,  t  e  A  , 

d  a  at 

are  as  defined  in  system  (5). 

The  mean  cumulative  delay  S  can  be  computed  from  the  solution  of  system  (8) , 


since,  by  definition. 


S  =  W, 


Before  proceeding  to  the  solution  of  systems  (5) ,  and  (8)  we  elaborate  on 
E  (u)d) .  With  reference  to  figure  3  we  define  the  following 

nd  :  the  number  of  packets  successfully  transmitted  during  the  first  CRI  of  a 

d-session  (Note  that  ndis  also  the  number  of  the  Poisson  arrivals  in  the 

enabled  arrival  interval  (t  ,  t 

n  n+l 

ud  :  the  sum  of  the  delays  of  the  packets  that  were  successfully  transmitted  during 
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k 


s 

t*. 


<r. 

►». 


r-' 


I 


I 


t: 


i 


the  first  CRI  of  a  d-session .measured  at  time  tn+^. 


v  :  the  sum  of  the  delays  of  the  packets  that  were  successfully  transmitted 
d 


during  the  first  CRI  of  a  d-session  after  the  first  CRI  begins. 
Using  the  above  defined  random  variables  we  can  express  as  follows 


“d  '  ud  +  vd 


if  1  <  d  <  A 


(10. a) 


u) .  *  u.  +  v.  +  n  (d-A)  if 
d  A  A  A 


d  >  A 


(10. b) 


Taking  expectations  in  (10),  and  using  known  properties  of  the  Poisson 
process  we  obtain 

,2 


Ad 


E  (o>d)  =  E  (vd)  +  2 


if  l<d<  A 


(11. a) 


AA 


E  (wd)  =  E  (v^)  +  +  X  A  (d-A)  if  d  >  A  (11.  b) 


Solving  systems  (5)  and  (8) 

Formally,  both  system  (5)  and  system  (8)  always  have  an  infinite  solution 


xd  =  +  d  e  A.  The  following  theorem  specifies  a  sufficient  condition  under 


which  both  systems  have  nonnegative  finite  solutions  (0<  x,  <  dcA  ,  d  <  +  <») 

—  a 

that  coincide  with  the  corresponding  algorithmic  sequences  induced  by  the  ZPDA, 
and  gives  upper  and  lower  bounds  on  such  solutions. 

Theorem  2 


If 


A  >  E  (Ifi) 


(12) 


then  (i)  system  (5)  has  a  solution  (yd)  d  e  A,  and  system  (8)  has  a 


solution  (z,}  d  e  A  ,  such  that 
a 


.* 


k  *>  V.  *> 

V.V.’W-'.A  v 


au  d  +  h  -  yd  -  au  d  +  6u  ,  d  e  A  (13) 

and 

y  d2+  <$  d+  C-  <  z .  <  y  d2  +  i  d  +  ?  ,  d  e  A  (14) 

'u  u  L  —  d  —  *u  u  * 

where  (*u,  3U»  3^*  Yu*  <SU»  Cu»  5^  are  real  coefficients  whose  expressions  are 
given  in  Appendix  A. 

(11)  the  algorithmic  sequences  {H^}  d  e  A  and  {W^}  d  e  A  coincide  with  the 

solutions  {yd)  d  e  A  and  {z^}  d  e  A  ,  respectively. 

The  proof  of  theorem  2  is  given  in  Appendix  A. 

Under  (12),  theorem  2  gives  the  following  bounds  on  H,,  W,,  d  e  A. 

a  a 

ad  +  et  t<\  £Hd  <  Hd(0)  4  aud  4  B„  ,  d  e  A  (15) 

Y  d2  +  i  d  +  ?,  <  W.  <  W*0)  Ayd2+6d+C  ,  d  e  A  (16) 

u  u  C  d  ~  d  —  d  'u  u  u 

The  above  ^initial  bounds  can  be  improved  using  the  method  described  below. 
Consider  the  case  where  we  want  to  improve  the  upper  bound  on  Hd.  To  this  end, 
define 

H^1}  £  E(£d)  +^2  pd(i)  h£0)-  h£0)+  E Ud)  +  au(EC£d)-d)  -(l+Xd)exp(-Xd)(au+fiu), 

l<d<A,  (17. a) 

Hd1>4  E(V  +  iIiPl1)  Hd-A+i  “  Hd°)+  EdA)-au(A-E(£A)),  d  >A,  (17.b) 

Under  (12) ,  the  coefficients  a  ,  3  have  been  chosen  so  that 


Given  a  natural  number  N  >2,  define,  for  n  >1, 


H(u)  A  E(^  +  ^  p^(i)  +  pd(i)  H(°)  f  l<d<A 

„(n) A  Ff/  x  .  £0  (i)  (n-1)  ®  (i)„(0)  . 

Hd  *  E(V  +  i=l  PA  Hd-A+i  +  ihn  PA  Hd-A+i  ’  A  d  N 


(19. a) 


(19. b) 


where  Nq  is  the  maximum  integer  satisfying  d-A+Ng<N 
For  n>l ,  combining  (17)  and  (19)  yields 

Hd0)  ==Hd1)  +  i^2  pd±)  (Hin_1)  ‘  Hi0)>  l£d<A 


(20. a) 


,(n)  _  „(D 


(n-1)  „(0) 


«d  '  =  +  ill  PA(i)  *  A<d<N 


(20. b) 


From  (17),  (18),  and  (19)  we  have  that,  for  every  d  e  fl,N)  the 

sequence  (Hy' }  is  a  monotonically  decreasing  sequence  of  upper  bounds  on 

d  iv>U 

H^.  The  upper  bound  ,  n>l,  d  £  [1,N)  can  be  computed  using  (20).  In  (20), 

A  is  taken  to  be  rational  so  that  d  assumes  finitely  many  values  in  the  interval 

[1 ,N) .  The  probabilities  p^(i)  ,  l£i<N,  l£d<A,  are  obtained,  in  a  recursive  form 

using  the  rules  of  the  algorithm  in  conjuction  with  straight  forward  combinatorial 

arguments;  their  expressions  can  be  found  in  [  16  }. 

To  improve  the  lower  boui.d  on  H^,  we  follow  the  same  procedure  except 

that  now  ^H,  is  used  instead  of  Hj^  ;  this  results  in  a  monotonically  increasing 
d  d 

sequence,  of  lower  bounds  on  H^,  for  d  e  (1,N).  The  same  method  can  be 

used  to  obtain  a  monotone  decreasing  sequence,  {Wjn^}n>Q  ,  of  upper  bounds,  and 
a  monotone  increasing  sequence  ,  of  lower  bounds  on  W^,  for  d  e  [1,N). 


/  / 


2.2.2  Mean  packet  delay  bounds  and  maximum  stable  throughput. 
From  theorems  1,  2  and  (1),  (9)  we  have  that,  under  (12), 


DQ  <  D  <  D 


where 


y  +6  +c»  Y  +5  +C 

Yu  *°u  Oil  -1  °  u 


n  A  1  _  1  _  Tu  u  „0  4  1  "l  u  u  -T« 

°  X  H^05  "  *(0tu  +Pu>  ’  X  X  (0tu  ■‘V 

The  bounds  given  by  (17)  can  be  improved  using  the  method  described  at  the 

end  of  the  previous  section.  The  improved  bounds  are 


where 


D  -  i 
n  X 


D  <  D  <  D 
n  —  — 


nn  1 _ ; 

*  X  (n). 


,  n  ■  l,2,i 


The  bounds  given  by  (21)  or  (22)  are  valid  for  all  X's  for  which  inequality 
(12)  is  satisfied.  In  Appendix  B,  it  is  shown  that  inequality  (12)  holds  if 

X  £  X  (m,A) 

where,  for  given  m,  X(m,A)  is  maximized  for  A  ■  A  .  The  values  of  X  (m,A  ) 

m  m 

* 

and  A  are  given  in  table  1,  for  different  values  of  the  parameter  m.  Table 
m 

2  gives  the  lower  bound  D^,  and  the  upper  bound  on  the  mean  packet  delay  D 

for  different  values  of  m  and  for  A  -  A  ,  where  A  is  a  rational  number  close 

m  m 

* 

to  A  ,  (selected  for  convenience  in  the  computation  of  the  bounds) . 


-  -  *  .  *■.  <V.*.  f.  .V.V.V.V .V.V.y  *.•  * IV 
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If  we  define  the  maximum  stable  throughput,  ,  of  the  ZPDA  as 
H  “  sup  {Xq  :  D  <  08  for  all  X  e  [0,Xq)}- 

then,  for  given  m,  X  (m,A*)  is  a  lower  bound  on  n>  since  (12)  is  only  sufficient 

m 

for  D  to  be  finite.  It  can  be  easily  shown,  however,  that  the  lower  bound  in 
(13),  increases  to  infinity  as  E(t^)  increases  to  A,  or  equivalently  as  X 

—  A  _  Ac 

increases  to  X(m,A  );  thus,  n*X(m,A  ).  A  more  intuitive  argument  can  also  be 

m  m 

used  to  show  that  (12)  is  a  necessary  condition  for  finite  mean  packet  delay. 

Since  the  algorithm  spends  on  the  average  E(A)  units  of  time  to  successfully  transmit 
the  packets  arrived  during  A  units  of  time,  inequality  (12)  is  necessary  for  the 
ZPDA  to  keep  up  with  the  arrivals. 

Finally,  note  that  to  account  for  the  overhead  created  by  the  m  control  bits 
per  slot  the  maximum  stable  thoughput  n  should  be  normalized  to 

_  n 

n  1  +  m  r/L 

where  L  is  the  length  of  a  standard  data  packet  in  bits,  and  r  is  the  length  of 
a  CMS  in  bits.  However,  since  the  feedback  from  a  CMS  is  of  the  E-NE  type,  r  can 
be  made  very  small — theoretically,  1  bit.  On  the  other  hand,  L  is  of  the  order  of 
1000  bits.  Thus,  the  reduction  in  throughput  due  to  the  overhead  is  small.  The 
same  normalization  should  be  applied  to  the  results  listed  In  table  2;  the 
values  of  X  should  be  divided  by  (1  +  m  r/L),  and  the  values  of  the  delay  bounds 
should  be  multiplied  by  (1  +  m  r/L) . 


3.  ADAPTATION  OF  THE  ZPDA  TO  THE  SATELLITE  CHANNEL 

Consider  now  a  channel  where  the  propagation  delay,  P,  is  much  larger  than 
the  packet  transmission  time.  A  typical  example  is  a  satellite  channel  where  the 
round-trip  propagation  is  equal  to  0.27  sec.  If  packets  of  1125  bits  are  transmitted 
with  a  rate  of  50  kbits/sec,  then  P  ^  12  slots.  Thus,  the  users  of  a  satellite 
channel  learn  the  outcome  of  a  transmission  P(»l)  slots  after  the  end  of  the 
transmission,  (i.e,  the  feedback  delay  is  equal  to  P+1  slots,  where  P>>1) . 

The  simplest  way  to  deal  with  the  problem  of  non-zero  propagation  delay  is  to 
treat  the  random-access  channel  as  P  =P+1  interleaved  zero-propagation  delay  channels. 
In  this  scenario,  slots  i,  i+P,  i+2P,...  of  the  actual  channel  form  slots  1,2,3,... 
of  the  i-th  interleaved  channel,  i=l,2,...,P.  The  Poisson  arrival  stream  is  sub¬ 
divided  into  P  interleaved  substreams  each  one  of  which  is  served  by  a  different  inter¬ 
leaved  channel  in  a  time  division  (TD)  manner.  Specifically,  a  new  packet  that  arrives 
during  slot  j-1,  j>l,  is  assigned  to  the  interleaved  channel  corresponding  to  slots 
j,j+P,  j+2P, . . . 

Suppose  now  that  the  ZPDA  is  executed  independently  on  each  of  the  P  interleaved 
channels.  Then,  clearly,  the  maximum  stable  throughput  of  the  channel  coincides  with 
that  induced  by  the  ZPDA,  i.e,  A(m,A*m).  The  induced  mean  packet  delay,  D^,  is  given 
by  (23), 

Dtd  =  P(D-  \)  +\  (23) 

where  D  is  the  mean  packet  delay  induced  by  the  ZPDA.  (23)  follows  from  the  fact 
that  on  the  average  a  new  packet  waits  half  a  slot  before  the  next  channel  slot 
of  the  interleaved  channel  it  is  assigned  to  begins,  and  P  (D-  -j)  channel  slots 
until  it  is  successfully  received  at  its  destination. 
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Bounds  on  are  obtained  by  substituting  the  bounds  on  D,  found  in  section 
2. 2. 2, in  (23).  Figure  4  shows  these  bounds  for  different  values  of  m,  and  for 
P  =  12  slots. 

Statistical  time  division  (STD)  is  an  alternative  rule  for  assigning  packets  to  the 
P  interleaved  channels.  This  rule  assigns  waiting  packets  to  the  first  available  inter¬ 
leaved  channel.  Specifically,  let  I  denote  the  instant  on  the  channel  axis  at  which 
a  CRI  ends  in  any  of  the  P  interleaved  channels;  let  also  s  be  the  time  on  the 
arrival  axis  for  which  it  is  known  that  all  packets  generated  in  the  interval 

[s,I  )  have  not  accessed  any  of  the  P  channels  yet.  Then,  all  packets  arrived  in 

s 

the  interval  [s,  s+min  (I  -s,A))  are  transmitted  in  slot  I  ,  i.e.,  they  are  assigned 

s  S 

to  the  channel  to  which  slot  Ig  belongs.  The  point  s  is  updated  as  in  the  ZPDA, 
i.e.,  at  each  Ig,  s-*-  s+min  (Is  -s,A);  initially  s=P,  since  at  the  beginning  of  the 
operation  of  the  system  and  for  the  first  P  slots  packets  are  assigned  to  the  P 

channels  according  to  the  TD  rule.  It  is  expected  that  the  STD  assignment  rule 
induces  improved  delay  characteristics  as  compared  to  those  induced  by  the  TD 
assignment  rule,  especially  when  the  traffic  is  light  to  moderate.  There  are 
two  reasons  that  support  this  claim.  First,  in  contrast  to  TD,  with  STD  no 
channel  is  idle  while  there  are  packets  waiting  for  their  first  transmission. 

Second,  since  the  arrival  axis  is  not  subdivided  a  priori,  as  in  the  TD  rule, 
the  CRl's  start  resolving  enabled  intervals  of  length  closer  to  the  optimal 
window  size  A.  However,  due  to  the  interdependence  in  the  operation  of  the 
interleaved  channels,  at  this  point  we  have  not  been  able  to  analyze  the  delay 
characteristics  induced  by  the  interleaved  system  operating  under  the  STD 
rule.  The  simulated  mean  packet  delay  performance  is  shown  in  figure  5  for 
m=16.  For  comparison  purposes,  in  the  same  figure  we  have  included  the  mean 
packet  delay  curve  for  ra“16  induced  by  the  TD  assignment  rule.  It  is  clear  that 


. V.  '---  -.V.- 


k  •  •  v  •  - 


_V  »  V. 


•V.V.W.V.V.V 


the  STD  scheme  offers  considerable  delay  performance  improvement  as  compared 
to  the  TD  scheme.  In  the  same  figure,  the  straight  line  at  25  slots  (two  times 
the  round  trip  propagation  delay  plus  one)  corresponds  to  the  minimum  mean  packet 
delay  induced  by  any  reservation  scheme.  As  it  can  be  seen  from  the  figure,  in 
the  satellite  environment  the  proposed  schemes  outperform  any  reservation 
scheme,  at  least  for  low  to  moderate  input  traffic  values. 
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Appendix  A 


Proof  of  Theorem  2 


Part  (i) 


Let  E  be  the  space  of  sequences  X  =  {x^}:  A+R.  Given  Y  eE,  define 


>’d3+1)-  bd  +  tfAc<it  ?t3)- deA-  J-0’1-2-- 


(A-l) 


where  b^,  c^,  d,teA,  are  defined  in  system  (5).  The  sequence  {Y^}, 
j=0,l,2,...,  generated  by  (A-l)  will  be  called  the  power  sequence  of  system 
(5)  with  initial  point  . 

Let  y^ ■  a  d+8  ,  deA.  To  establish  the  existence  of  a  solution 
d  u  u 

Y={y .} ,  such  that  0<y,  <a  d  +  8  ,  deA  it  suffices  to  choose  a  ,8  so  that 

J  A  *  ~^d  —  U  U  >■  ii 


o<y^^  f°r  every  deA 


(A-2) 


Under  (A-2)  and  since  b  j>  0,  c  _>  d,  teA,  we  have,  that 


0  <  y^+1)  <  y^,  deA,  j=0,l,2,... 

—  n  —  H 


(A-3) 


Thus,  under  (A-3),  the  following  limit  exists 


lim  y^  =  y,  ,  deA 

j-*»  a  a 


(A-4) 


Passing  to  the  limit  as  j-*»  in  (A-l)  shows  that  the  sequence  Y  solves 


system  (5) 


From  (A-3) ,  (A-4)  we  have  that 


0<y,<ad+8  ,  deA 

—  d  —  u  u 


Next  we  choose  a  ,  8  so  that  (A-2)  is  satisfied, 
u  u 
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After  straightforward  manipulations,  by  (A-l) ,  we  obtain 

yd1>=  yd°)+  E(V  +  “u  (E(£d)  _d  -Pd(1))"  6u  Pd(1)  »  1  1  d  1  A  (A-5.a) 

yd1)=  yd°>+  E(V  '  au  (  A  _  E('eA))  »  d  >  A  (A-5.b) 


From  (A-5.b)  we  see  that  (A-2)  is  satisfied  only  if  condition  (12)  in  the 

theorem  holds.  Under  (12), it  can  be  readily  seen  from  (A-5.a)  that  (A-2)  is 

satisfied  if  we  choose  a  ,  8  as  follows 

u  u 

E(V 

°u  "  7uI77  ■)  »  "  8UP  (♦(d)) 

u  A- E(£a)  u  1<d<A 


where  \p  (d) 


E(£d}  +  “u  (E(£d}  “  d  '  Pd^1)) 

i^i) 


,  and  where  Pd(l)  =  (1+Xd)  exp  (-Xd) . 


construct  a  lower  bound  on  yd  we  let  yd  =*  ot^d  +  6^,  deA,  and  we  choose 


and  6^  such  that 


,<°><  y<»  ,  deA 


(A-7) 


From  (A-5)  it  can  be  readily  seen  that  (A-7)  holds  if 


ot»  -  a  ,  Bp  “  inf  (<|>(d)) 
l  u  Kd<A 


where  au>  4>(d)  are  as  given  above. 

The  expectations  E  (£d) ,  1  £  d  £  A,  needed  for  the  computation  of  the 
coefficients  of  the  linear  bounds  can  be  computed  as  shown  in  Appendix  B. 


System  (8) ;  similarly  to  system  (5),  the  solution  Z  -  {zd)deA  to  system  (8) 
may  be  obtained  as  the  pointwise  limit  of  its  power  sequence  {Z^},  j-0,1,2,.... 


defined  by 


l- 

k 


v."  .1.  ....  ,T  .-.  1  ..»■  T  ■ ^  ^  w.  • 
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,(j+D 


=  b'  +  £  c  z 

d  teA  dt  d 


(j) 


,  deA 


(A-8) 


provided  that  z^  is  such  that 
d 


>  0,  z<1}  £  z<°>  ,  for  every  deA  (A-9) 


Under  ( 12)  ,it  can  be  shown  by  direct  substitution  in  (A-8),  that  (A-9)  holds  if 


we  choose  z^  =  Y  d2  +  <$  d  +  C  ,  where 
d  u  u  u 


Y„  = 


AA 


u  2(A-E(£| A) ) 


5  = 


E(o)a)  -\A2  +  Yu  (A2  +  E(£2)  -2  A  E(£|A)) 


A  -  E(£ | A) 


C  =  sup  (<}>(d) ) 
l<d<A 


4»(d)  = 


E(u  )  +  Y  .(E  (lh  -d2)  +  <5  (E(£  )  -d) 
dud  u  a 


(1+Xd)  exp(-Xd) 


(Y  +  A  ) 
u  u 


Under  (12), and  for  every  deA,  the  sequence  { ^ } j>q  is  non- increasing;  thus 


z  <  y  d  +6  d+£  . 
d  —  u  u  u 


The  construction  of  a  lower  bound  Y^d  -f-d^d  on  z^  is  similar  to  that 
of  the  upper  bound;  the  coefficients  Y^»  and  are  given  below 


Y/  =  Y„’  *  5.,’  *>9  =  inf  (♦(<*)> 

U  c  U  C  l<d<A 


The  expectations  in  the  expressions  of  the  coefficients  Y  .  A  »  C  *  can  be 

u  u  u  c 

computed  as  shown  in  Appendix  B. 

Part  (ii) 

The  proof  involves  two  steps.  In  step  1  we  show  that  the  solutions  found  in 
part  (i)  are  unique  in  the  class  of  sequences 


\  ■  . 


\  ».  \  \  \  \  •-  *v  •*, 


In  step  2,  we  show  that  the  algorithmic  sequences  {H^}  deA,  and  {W^}  deA 
belong  to  E^,  and,  therefore,  coincide  with  the  corresponding  unique  solutions 
in  E2- 

Step  1:  As  in[17],  we  shall  call  the  system 


-  Bd +  4catV deA 


majorant  of  the  system 


(A-10) 


x 


d 


bd  + 


1  c, 

t£A  dt 


V 


deA 


(A-ll) 


if  the  following  inequalities  hold: 

|cdt |  -  Cdt’  bd  -  Bd  *  d»teA 

A  solution  of  a  system  of  type  (A-ll)  that  is  a  pointwise  limit  of  its  power 
sequence  with  initial'  point  -  o,  is  called  its  principal  solution. 

With  the  above  definitions,  we  state  Lemmas  A.l,  A. 2  below,  which  are  essentially 
theorems  I,  II,  &  2  of  [17]. 

Lemma  A.l  If  the  majorant  system  (A-10)  has  a  nonnegative  solution  S,  then  both 
system  (A-10)  and. (A-ll)  have  principal  solutions  S*,  S*,  respectively.  Moreover, 

0  <  l8$l  1  1  8d  . 

Lemma  A. 2  If  the  majorant  system  (A-10)  has  a  nonnegative  solution  S,  then  the 
principal  solution  S*  of  system  (A-ll)  is  unique  in  the  class  E(S*)  C  E,  defined  as 


follows 


E(S*)  =  {  X  e  E  :  sup 


deA  s-  * 

d 


Furthermore,  S*  is  the  pointwise  limit  of  any  power  sequence  of  (A-ll) ,  with 
initial  point  any  point  in  E(S*) 

Using  the  above  lemmas,  we  shall  show  that  both  the  solution  Y  ={y<j}  of 
system  (5),  and  the  solution  Z  *{  jjeA  sYstem  (8)  are  unique  in  £2*  We  start 
with  system  (8).  Since  system  (8) has  a  nonnegative  solution,  Z,  and  is  majorant  of 
itself,  from  lemma  A.l,  it  has  a  principal  solution  Z*  ■  {z*}  deA  *  According  to 
lemma  A. 2,  the  solution  Z*  is  unique  in  the  class 

l*dl 

E(Z*)  *  {XeE  :  sup  — 5—  «»} 
deA  2d 

provided  that  Z^e  E  (Z*) .  Since  by  definition  (see  part  (i)),  Z^e  E2,  the 
solution  Z  will  be  unique  in  E2>  if  E(Z*)  =  E2*  To  show  that  the  latter  holds,  we  use 
the  following  lemma  whose  easy  proof  is  omitted. 

Lemma  A. 3  Let  F  »{fd)d£.^,  G  two  sequences  in  E.  If 


d  deA 


(a)  fd  >0,gd  >0,  deA 

fd 

(b)  sup  —  <°° 

deA  gd 

(c)  inf  —  >0 
deA  8d 


sup  —  ■  <00 

deA  8d 


iff  sup 


deA  £d 


,  X  e  E; 


lxd' 

i.e. ,  the  classes  E_  *  (X  e  E:  sup  — - —  <°°}  and  E„  ={X  £  E: 


>»Vv 


Let  {Z  },  n=0,l,2,...,  be  the  power  sequence  of  system  (8),  with 
z^^=0,  deA.  Clearly 

bd  ^  e  >  0,  for  every  deA.  (A-12) 

Also,  it  can  be  shown  by  induction  that, 

z<n)  =  n  E(u>d)  -  n  (n-1)  X  A(A-E(^A))  (A-13) 

for  every  deA,  n  _>  1,  such  that  d  >  n  A.  For  d  >  2  A,  let  n  =  |^j-  1  in 
(A-13),  where  [xj  denotes  the  maximum  integer  not  exceeding  x.  Then,  using 
the  inequalities  E(wd)  >  X  A  (d-A)  +  ^~  ,  (see  (11. b)),  L  ^ -1  >  ^  and 
A  L  ^  J  1^1®  (A-13)  yields 

S«>  |  x  d2  -  X  A  d  ,  d  >  2  A  (A-14) 

From  Lemma  A.l  and  the  fact  that 


b'A  _>  0  ,  cdt  _>  0  ,  d,  teA,  we  have  that 


z*  -  ^d11^  -  0  »  deA»  n  -  1 


From  (A-12)  and  (A-15)  we  obtain 


Z*  -  e  >  0  »  d£^ 


From  (A-14),  (A-15),  and  (A-16)  we  conclude  that. 


z  ^ 

inf  d 


deA  .2 


>  0 


From  (14), and  the  inequality  z*  £  z^,  (lemmaA.l),  we  have 


Z," 

sup  d 

deA  d2 


(A-15) 


(A-16) 


(A-17) 


(A-18) 


Finally,  from  (A-17),  (A-18),  and  lemma  A. 3  we  have  that  E(Z*)  =  • 

To  establish  the  uniqueness  in  of  the  solution  Y  -  (yd)  d£^  of  system  (5) 


' *.  .v.v  V 
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we  proceed  as  follows.  Comparing  Che  constant  terms  of  system  (5)  to  the 
constant  terms  of  system  (8)  we  can  readily  show  that  there  exists  a  positive 
number  M<°°,  such  that 

£  Mb^  ,  for  every  deA  (A-19) 

From  (A-19)  we  have  that  system  (A-20) 

xd  ■  “d  +  txAcdt  xt  •  d£A  <A-20> 

is  majorant  of  system  (5).  Also,  note  that  the  sequence  Z  *  <V  deA’  *d  '  Mzd’  deA’ 

solves  system  (A-20),  since  Z  ■  {z.}  ’  is  a  solution  of  system  (8).  Furthermore, 

a  deA 

system  (A-20)  has  a  principal  solution  1  «  Cz*}  ^  with 

-  M22  ,  deA  (A-21) 

where  Z*  =  {z,}  is  the  principal  solution  of  system  (8) 
a 

Since  system  (A-20)  is  majorant  of  system  (5),  from  lemma  A. 2,  we  have  that  the 
solution  Y  ■  (yd)  is  unique  in  the  class 

E(Z*)  *  {XeE  :  sup  -Lill-  >  ®  } 
deA  x* 

d 

provided  that  e  E(Z*) .  From  (A-21)  we  have  that  E(Z*)=E(Z*).  However,  it  has 

already  been  shown  that  E(Z*)  -  E2 ;  thus,  E(Z*)  =  £2*  Since,  by  definition  (see  part 
(i)),  Y^  eE2t  we  conclude  that  Y  is  the  unique  solution  of  system  (5)  in  £2* 

Step  2 

Here,  we  show  that  the  mean  session  length  sequence  {H^}  deA  belongs  to  the  class  E2. 

Define  the  random  variables  h ?  •  min  (h.,n),  n-1,2,3,...,  deA,  where  h,  is  the  length 

d  d  c 

of  a  d-seaslon.  Clearly, 


I 


fl< g* )  deA 


(A-22) 
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If  we  let  H*  =  E(h j} ,  then  from  the  definition  of  hj  and  the  operation  of  the 
d  d  d 

algorithm  we  have  that 

Hd  i  bd  +  teA  cdt  ■  deA  (A-23) 

Also,  since  h*j  £  n  a.e.,  deA,  the  sequence  (H^)  d£^  belongs  to  E^*  f°r  every  n. 
Using  (A-23)  and  the  fact  that  Y  -  {yd>  d£/^  is  a  solution  of  system  (5)  that  belongs 
to  it  can  easily  be  shown  that 

H"  <  y .  ,  deA,  n-0,1,2,...  (A-24) 

a  —  a 

From  (A-22)  and  the  monotone  convergence  theorem  we  have 

Hj /  H,  ,  deA  (A-25) 

a  a 

From  A-24)  and  (A-25)  we  conclude  that  Hj  <  yd;  thus  (Hd)  de^  e  E^.  Since 
{H^}  d£^  solves  system  (5)  it  coincides  with  its  unique  solution  Y  in  E^. 

The  proof  for  the  mean  delay  sequence  {Wj}  deA  is  similar. 
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Appendix  B 


Computation  of  E(£_,)  ,  E (■£ ,)  and  E(v  ,) 


Let  X,  denote  any  of  the  random  variables  t , ,  or  v,,  and  define 
a  add 

E(Xj|k):  The  conditional  expectation  of  the  random  variable  X^,  given  that  the 
enabled  arrival  interval  contains  k  packets,  k  ^  0. 

The  conditional  expectation  E(X,  |k)  does  not  depend  on  d;  therefore,  we  may 


write 


E(xd)  =  kl0  E<x1|k)  e 


-Xd  (Xd) 
e  k! 


(B.l) 


Note  that,  since  the  length  of  an  enabled  interval  is  at  most  A,  we  have  that 

E(XJ  =  E(X.)  ,  for  d  >  A 
d  A 

Define  Lfc  =  ECCjk),  ^  -  ECeJjk),  and  =  E(vjk),  k  >  0.  Using  the  rules  of 
the  algorithm,  the  quantities  LR,  Mj,  and  VR,  k  _>  0,  can  be  computed  recursively. 


as  follows. 


where 


M0  -  Mi  -1 


LQ  -  Lx  -1  ;  Lk  -  1  +  m  ^  b^  (k,l/m)  L ,  k  >  2, 


b  (k,l/m)  -  <k)  (±)J  (1-  i)k  j 
j  j  m  m 


k-1  k-k. 


Mk*l4m  jS^bj  (k,l/m)Mj  +  2(Lk-l)  +  m(m-l)  ^  ^  q^kj.kj,)  L^L^,  k  >  2, 


where 


qk(kl,k2)  “  k1!k2!(k-k1-k2)!  (m 


k.  +k  k-k 

r  (m}  1  2  (1"  b  1 


vo  ■  °*  V1  ■  1 


\  ’  k+1L  bj<k'1/n)  vj  +  Kl0  j!0  i£i  bj(k'I>  bi(k-J'5ic>  Li  - 


k  >  2 


Using  the  above  recursions,  a  finite  number,  M,  of  terms  from  the  infinite 
series  (B.l)  can  be  easily  computed.  Also,  for  large  K  values,  and  based  on  the 
recursive  expressions,  simple  upper  and  lower  bounds  on  E{X^|k}  can  be  developed 
Those  bounds  can  be  used  to  tightly  bound  the  sum 


00 
£ 

k-Mfl 


(Xjk) 


-Xd  <Xd)‘ 
®  k! 


The  condition  A  >  E(-£. ) . 


Given  the  algorithmic  parameters  A  and  m,  and  in  view  of  (B.l),  the  condition 

A  >  E(£^)  in  theorem  2  can  be  equivalently  expressed  as  X  <  X(m,A) ,  where 

T  -XA  (XA)k 


00 


} 


m 

A(m,A  ) 
m 

* 

A 

m 

2 

0.5021 

2.739 

3 

0.5782 

3.182 

4 

0.6168 

3.665 

5 

0.6424 

4.347 

6 

0.6621 

5.422 

7 

0.6785 

6.476 

8 

0.6924 

7.179 

10 

0.7145 

8.051 

12 

0.7317 

8.636 

14 

0.7456 

9.091 

16 

0.7572 

9.492 

18 

0.7672 

9.847 

20 

0.7758 

10.174 

22 

0.7834 

10.479 

24 

0.7902 

10.766 

26 

0.7963 

11.025 

28 

0.8019 

11.303 

30 

0.8070 

11.555 

32 

0.8116 

11.801 

Table  1  The  maximum  stable  throughput 
_  * 

A(m,A  ),  and  the  optimal  window  length 
*  m 

A  ,  for  the  ZPDA. 
m 


m  "  8 


m  *  16 


0.1 

1.7068 

1.7068 

1.6820 

1.6820 

1.6724 

1.6724 

0.2 

1.9945 

1.9945 

1.9233 

1.9233 

1.8958 

1.8958 

0.3 

2.4494 

2.4494 

2.2749 

2.2749 

2 . 2088 

2 . 2088 

0.4 

3.3041 

3.3041 

2.8426 

2.8426 

2.6795 

2.6795 

0.5 

5.5664 

5.5665 

3.9299 

3.9299 

3.4589 

3.4589 

0.6 

33.8550 

34.1730 

7.1034 

7.1040 

5.0305 

5.0305 

0.65 

- 

- 

- 

- 

6.7409 

6.7414 

0.68 

- 

- 

44.1240 

44.6640 

- 

- 

0.7 

- 

- 

- 

- 

11.1590 

11.1870 

0.74 

- 

- 

- 

- 

32.1660 

33.0430 

Table  2  The  lower  bound  and  the  upper  bound  D  on  D, 
for  the  ZPDA. 
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Figure  4 

The  mean  packet  delay  D^, 


for  P  -  12 
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